Steady state entanglement in the mechanical vibrations of two dielectric membranes 
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We consider two dielectric membranes suspended inside a Fabry- Perot-cavity, which are cooled to 
a steady state via a drive by suitable classical lasers. We show that the vibrations of the membranes 
can be entangled in this steady state. They thus form two mechanical, macroscopic degrees of 
freedom that share steady state entanglement. 
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Introduction - Optomechanical systems in which 
electromagnetic degrees of freedom couple to the mechan- 
ical motion of mesoscopic or even macroscopic objects 
are promising candidates for studying the transition of a 
macroscopic degree of freedom from the classical to the 
quantum regime. These systems can also be of consid- 
erable technological use, e.g. for improved displacement 
measurements [l| and their application in the detection 
of gravitational waves 0] ■ 

Optomechanical devices have therefore attracted con- 
siderable attention in recent years and micromirrors have 
been cooled by radiation-pressure In many setups, 
one of the end mirrors of a Fabry-Perot cavity undergoes 
a mechanical vibration and the coupling between cav- 
ity photons and the mirror motion emerges because the 
resonance frequency of the cavity depends on its length 
and hence on the position of the mirror. Recently devices 
have been introduced in which the motion of a membrane 
that is inserted into a Fabry-Pcrot cavity formed by rigid 
mirrors couples to the cavity mode 0, . 

Whereas the ground state and hence the quantum 
regime has not yet been reached in experiments, this 
has been predicted to be achievable if the mechanical 
oscillation frequency is larger than the cavity linewidth 
[1], a regime that has recently been observed In the 
quantum regime, it is then interesting to explore entan- 
glement in mechanical i.e. macroscopic degrees of free- 
dom [1, [^, Possibilities to entangle the motion of 
a cavity micromirror with the electromagnetic field in 
the cavity have thus been explored in various approaches 
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FIG. 1: The setup: Two mechanically vibrating membranes 
(brown) are suspended inside a Fabry-Perot-cavity, which is 
driven by external lasers (green). Dissipation occurs via me- 
chanical damping (brown) and cavity decay (blue). 



to very low temperatures [13[. In contrast, our scheme 
generates steady state entanglement between mechanical 
degrees of freedom by cooling them via radiation pressure 
which only uses classical light sources. Our approach is 
not restricted to the specific setup mentioned here but 
also applies to other devices with optomechanical cou- 
plings [3(1, between two mechanical and two cavity modes. 

Model - We consider a Fabry-Perot cavity with two 
dielectric membranes in its interior (c.f. figure [1]). In 
this setup the optical resonance frequencies of the cav- 
ity depend on the positions of the membranes and can 
be derived from the boundary conditions of the field 
in the cavity 17|. Let — 3L (3L) and qi ((72) be the 
positions of the left (right) rigid mirror and the left 



Here, we consider a Fabry-Perot cavity with two di- 
electric membranes suspended in its interior (c.f. figure 
[T|) and assume that two cavity resonances arc driven by 
external lasers. With suitable lasers the mechanical vi- 
brations of the membranes are cooled and asymptotically 
driven into a steady state. We show that the mechani- 
cal vibrations of the two membranes can be entangled in 
this asymptotic state. Entanglement between mechanical 
oscillators has been discussed previously but was either 
only found in the transient regime [lH (not the steady 
state) or required either a drive with non-classical light 
[1^ or mechanical oscillators that had been pre-cooled 



(right) membrane and T the transmissivity of the mem- 
branes. We denote the field modes in the left, center 
and right part of the cavity by ui, U2 and U3. For a 
mode with wavenumber fc, the boundary conditions read 

= 0, ui(gi - 0) = U2{qi + 

g-j rfMi(<;i-0) _ du2(gi+0) _ 



U3(3L) 
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dU3{q2 



dz 



dz 



PJ, ui{-3L) 
0), ^2(92 - 0) 
fc?/m(gi) and 

r] = 2^(1 - T)/T. In analogy to 
ansatz ui = Ashi{k{q + 3L)), U2 = B cos{kq) 
and U3 = C sin(fc((j — 3L)) and obtain the transcen- 
dental equation, [cos(3fcL) — 77 cos(fcgi) sin(fcgi -|- 3kL)] 
[sin(3fci) + 1] sm{kq2) sm{kq2 — 3kL)] + [cos(3fcL) -|- 



= kriU2{q2), where 
0, we make the 
B sm{kq) 
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T]cos{kq2) sin(fcg2 — 3kL)] [sin{3kL)+r]sm{kqi) sin{kqi + 
3kL)] = 0, from which the optical resonance frequencies 
bj = kc [c is the speed of hght) can be found. 

To obtain optomechancial coupUng between two op- 
tical and two independent mechanical modes (see be- 
low), we choose the equilibrium positions of the mem- 
branes, (7oi and (702, to be goi = and go2 = 2L. For 
non- vibrating membranes, the optical resonance frequen- 
cies of the membrane cavity systems are then given by 

, , n-KC , , wnc 9c , c ^^^ — 1 i cosO\ , , 

^an - — , (^bn - — 2L ~^ 2L '^'^^ [ 2~ ) ' ^b'n — 

ni^c _ Be c_ ( cos e N j _ m-KC ttc 0c 

L 2L 2L '-^'^ \ 2 ; ' '^''^^ "^cm — 3L 6L 3L ' 

where n and m are positive integers. For membranes with 
low transmissivity, the frequencies ujan , ^bn and Wcm with 
m = 3n lie close together whereas ujb'n is separated from 
this triplet and we thus focus on Wan, i^bn and ujcm- 

In the case of vibrating membranes, the optical reso- 
nances depend on the motion of the membranes and their 
frequencies become functions of gi and 52, e.g. ujaniqi, Q2) 
(The assumption that the optical resonances only depend 
on the membrane positions, not their momenta, is only 
valid if the membrane oscillations are much slower than 
the optical round-trip time, i.e. ujm ^ \ujx — LUy\ for 
x,y = an,bn,cm which we confirm below.). To ob- 
tain these functions, we write them as a power series up 
to linear order in the membrane positions, LL!x{qi,q2) = 
i^xiqoi, qo2)+^xi{qi-qoi)+^x2{q2-qo2) for X = an,bn,cm 
iloi = —L and qo2 — 2L), expand the transcendental 
equation up to linear order in qj — qoj {j ~ 1,2) and 
solve it for zeroth and linear order separately to obtain 
<^a;(goi, 902) = i^x,^xi and ^2,2 for x = an,bn,cm. For 
our choice of the membrane rest positions, qoi — —L and 
go2 = 2L, the mode Wan does not couple to the membrane 
motions, ^ai = Ca2 = 0, and we discard it. The other cou- 
phngs are 61 ^ + ^,T) ^, 62 « (| - §,T) ^ 
and = -^,2 « - - WE^) for T « 1 and 
n, m 3> 1. Higher order terms in the expansion of the 
frequencies give rise to additional coupling terms, also 
for Uam but these are negligible compared to the lin- 
ear couplings. For symmetric membrane rest positions 
(goi = -i, qm = L), one would get Cbi+62 Cci+^c2 = 
and the photons would only couple to the breathing 
mode, qi — q2, whereas the center of mass mode, qi +q2, 
would not be cooled. Note also the two optical modes 
are needed to cool two mechanical modes. 

The corresponding Hamiltonian that describes the mo- 
tion of the membranes and the cavity modes reads 

j — 1,2 x—bn,cm 

+ (^x+ ^-iii) (1) 

x—bn,cm j — 1,2 

where pj and qj are the momentum and position of mem- 
brane j (j = 1,2). Both membranes have the same effec- 
tive mass m and mechanical resonance frequency and 



the optical modes with creation(annihilation) operators 
o\jn{^bn) and aj„j(acm) are driven by classical lasers with 
Rabi frequencies rif,Ti and ^cm- We have redefined the po- 
sition variables qj — q^j — > qj and write the optical modes 
in frames that rotate at the frequencies of their respective 
driving lasers, A^, = lOxo - w^c.Laser [11] ■ In eq. H]), we 
have also assumed that each laser only drives one cavity 
mode, which sets an upper bound to the permissible Rabi 
frequencies, |$lb„|, |17cm| < \uJbn -Wcm| ~ j^c\/T / L (to 
leading order in T <C 1). This in turn limits the amount 
of entanglement that can be generated. 

The linear optomechanical couplings, i,xjqj can be ex- 
ploited to cool the membranes and drive them into a 
steady state. They will furthermore generate entangle- 
ment between the mechanical vibrations via the optical 
modes as we will show. 

Equations of motion - Taking into account cavity 
decay and mechanical damping of the membranes, the 
Hamiltonian ([1]) gives rise to the Langevin equations [l^ , 

q-j ^UmP] , P'j = -t^m^j - ij^Pj ~2^(,xjala^ + Q, (2) 

X 

where dots denote time derivatives and [•]* a complex 
conjugate, a™ and Q are the optical and mechanical 
input noises and Tx and 7 cavity decay and mechan- 
ical damping rates. The relevant nonzero correlation 
functions of the noise operators are {a]^{t) (a™)^ {t')) = 
SxyS{t — t') for x,y ~ bn,cm and (Cj(^)Ci(^')) = 
^ {2n^^^ + I)dji6{t — t') for j,l = 1,2, where n^^ = 

(exp(?ki;,„/A;BT) — 1) ^ is the thermal phonon number 
of the mechanical environment at temperature T, ks is 
Boltzmann's constant and we have assumed fc^T ^ fiWm- 
We split the operators in ^ into their steady state 
expectation values and quantum fluctuations, ~ Cx + 
5x , Pj = Pj + 5pj and qj = Qj + Sqj . The constant steady 
state expectation values are given by the equations 

where i^x = '^x + CxiQi + £,x2Q2- 

We are interested in a regime of high photon num- 
bers in the cavity, in which the steady state expecta- 
tion values are much larger than the quantum fluctua- 
tions. In this regime we can neglect all terms of higher 
than linear order in the fluctuations, 6x, Spj and Sqj 
in ^ (We have confirmed this approximation numeri- 
cally.). The asymptotic state of the quantum fluctua- 
tions for the linearized equations is then a zero mean 
Gaussian state which is fully characterized by its covari- 
ance matrix Vij = 2Re((Oi — (Oi)) {Oj — (Oj))), where 
O = iSqi,Spi,Sq2,dp2,Xbn,Yb„,Xcm,Ycm) with Xx = 
(Sx + Sl)/V2 and Yx = -i^Sx - Sl)/V2. We solve the 



3 




FIG. 2: The entanglement of the two mechanical vibrations 
in the steady state as measured by the logarithmic negativity 
as a function of Cbn and Ccm. The phases of Cbn and Ccm 
do not affect the entanglement. Left plot: Atn = 4.07MHz, 
Acm = 20.84MHz and T = 8mK = 1000). Right plot: 

At„ = 6.12MHz, = 33.18MHz and T = lOOmK (n^,„ = 
13085). The remaining parameters are Um ~ IMHz, m = 
10"^g, T = 0.2, L = 1mm, n = 2 x 10^ = 3n = 6 x 10^ 
Qoi = -L, qo2 = 2L, Fbn = Tcm = t-^m/lO, and Q = 10^ 
for both plots. A(,„, Acm, Ci,„ and Ccm have been optimized 
numerically for each case. 



linearized Langevin equations for the fluctuations to ob- 
tain the steady state covariance matrix V in the same 
way as in From V, the steady state entanglement 

as measured by the logarithmic negativity En {En 
means the state is entangled) can then be computed [l6| . 

Steady state entanglement - We consider an example 
where both membranes have a transmissivity T = 0.2, 
an efi^ective mass of m = lO^^g and a mechanical res- 
onance frequency of ujm. = lO^Hz The cavity is 
6mm long, hence L — 1mm. For driving lasers of about 
lOOOnm wavelength, the closest cavity modes have num- 
bers n = 2 X 10'^ and m = 3n ~ 6 x 10"^. For these pa- 
rameters, the optomechanical couplings attain the values 
6n,i = 1.90kHz, ^b„^2 = 6.75kHz and ^cm.i = -Lm.2 = 
—4.53kHz (Note that we work in units, where Sqi and dq2 
are dimensionless and given in multiples of y^h/{muj„i)). 

Cooling to the quantum mechanical regime is pos- 
sible if the mechanical oscillation frequency is larger 
than the optical lincwidth [^, 0| and we thus assume 
Fbn = = ujm/W. Thc mcchauical Q is taken to be 
Q = 10^, consistent with 0]. For thc mechanical envi- 
ronment, we consider two temperature values, T = 8mK 
(n^^ = 1000) and T = lOOmK {n^^ = 13085). 

Figure [2] shows the entanglement of the two mechani- 
cal vibrations in the steady state measured by the loga- 
rithmic negativity, E^, as a function of the steady state 
electromagnetic fields in the cavity, ctn and Ccm- Since 
the Hamiltonian ([1]) only contains the photon numbers, 
the entanglement is insensitive to the phases of fields in 
the cavity and thus also to the phases of Cb„ and Ccm- 
The linearization of eq. ([2]) requires |c6ra|,|ccm| 3> 1. 
In the left plot, we have |cf,Ti| ~ 30, but we also find 
En = 0.195 for \cbn\ = 60, |c„„| = 386.4, Abn = 4.2MHz 
and Acm ~ 20.9MHz. The values for Q and T, we as- 
sume here, are currently hard to achieve simultaneously 
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FIG. 3: The steady state entanglement of the two mechanical 
vibrations for T = lOOmK, c^n = 60 and Ccm = 486. a: En 
as a function of A^n and Acm, b: En as a function of Tt,^ 
and Fcm, c: En as a function of m and cjm and d: En as a 
function of Q and T. All other parameters are the same as in 
figure [2l 

but the entanglement persists in a larger parameter range 
as shown in figure [3ji. 

The Rabi frequencies Qbn and flcm that are needed to 
generate the values of Cf,„ and Ccm in figure[l]are less than 
llGHz (left plot) and 23GHz (right plot). The diff'erence 
between the resonance frequencies LUbn and Ucm on the 
other hand is \ujbn ~ujcm\ ~ 57GHz and the separation of 
these two modes from other resonances is much larger, so 
that thc lasers indeed only drive one resonance mode as 
assumed in eq. ([1]). |a;h„ — Wcml oc L^^ so that |a;h„— Wcm| 
would even be larger for a shorter cavity. By reducing 
the cavity length one could thus employ stronger driv- 
ing lasers and create substancial entanglement even at 
higher T. The input laser powers Px are related to the 
Rabi frequencies by = fiWa;,Laser|f^Kp/(4r^), which 
implies that laser powers between 0.6 pW and 60 /iW 
arc required. Furthermore tUm *C \uJbn — ^cm\ and the 
derivation of the cavity resonances and consequently the 
form of Hamiltonian ([1]) arc well justified. 

Thc steady states are furthermore characterized by 
phonon numbers of the vibration fluctuations, nj = 
^{Sp^^ +5q^ ^l), of ni < 3,n2 < 5 for T = 8mK 
and rii < 5,n2 < 10 for T = lOOmK, and an entropy 
of thc reduced density matrix of the vibrations, Sm = 
-Trphotons(plog2/c) of less than 4, < 4, for both, 
T ~ 8mK and T = lOOmK. The steady state is thus 
indeed in thc quantum regime. The linearized Langevin 
equations can be cast in matrix form O = OA + n, where 
n is the vector of the noises [l3| • We have conflrmed that 
all eigenvalues of A have negative real parts below IkHz 
which ensures that there is a unique steady state that is 
reached within milliseconds. 

To further corroborate the robustness of the entangle- 
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ment, we studied its dependence on fluctuations in the 
driving lasers and on variations in several system param- 
eters. The results for T = lOOmK are shown in fig. [31 
Fig. [3^ shows the dependence of Ejy on the detunings of 
driving lasers, A^n and Acm (c^n = 60 and Ccm = 486), 
fig. [5)3 the dependence on the cavity decay rates r;,„ and 
Tcm, fig- [3fc the dependence on the membranes resonance 
frequency, w,„ and its effective mass, m and fig. [Sjl the 
dependence on the mechanical Q and the environmen 
temperature T. In all cases there is a substantial pa- 
rameter region which shows entanglement. Note that we 
assign En = to all points where there is no well defined 
steady state due to heating, i.e. where an eigenvalue of A 
has a positive real part. Furthermore, our driving fields 
are optimized for the values in figure [2] For different val- 
ues of Tbn, Fcm, Q, 7, Or m, slightly modified laser 
drives will yield more entanglement. 

One can obtain some intuitive indications why the 
membrane vibrations become entangled. The Hamil- 
tonian corresponding to the linearized Langevin equa- 
tions for the fluctuations Spi , 6qi , dp2 , Sq2 , Sbn and 

is H ^ i^j2j{Spj + Sq]) + E.i^As. + 

X i^xjCx^qj^i + h.c). In the parameter regime of in- 
terest, we have fj,x ^ ICaiiCxI and the photon degrees of 
freedom can be adiabatically eliminated to obtain the 
effective Hamiltonian H = ^{5pl + 5pl) + ^^^^^Sql + 
i^l^Sq^+!^SqiSq2, where J/, = -2 Cxj|c.| V/^. and 
1^12 -'iJ2x^xi^x2\cx\'^ / ^J■x■ The ground state of H and 
hence states close to it are entangled in regimes where 
1 1^12 1 ~ oJm or larger, which is the case for the parameters 
in figure [5] 

Entanglement verification ~ To verify the created 
entanglement quantitatively in an experiment, several 
quadrature correlations need to be measured [20| . 
This may be achieved by employing at least two fur- 
ther weak probe lasers, similar in spirit to the scheme 
in which drive cavity modes that do not partici- 
pate in the entanglement generation, i.e. modes with 
71 7^ 2 X 10^ or m ^ 6 X 10"^. The equation of mo- 
tion for the fluctuations of the probe field, 6y, reads 

4 = -itJ'ySy + a/T^'^" ~ + ^1/2(^2) in a 

frame that rotates at the frequency of the probe laser, 
lol- Here Sqj ~ "^(Qf + ^j) ^nd we have assumed 
\S,jCy\ <^ ujm and applied a rotating wave approxima- 
tion. In a Fourier transformed picture in the labo- 
ratory frame this equation reads —iiuj — ^y)5y{Lj) — 
yr;4"(a;) - i^[^y^Ci{oj - ujl) + Cy2C2iuj - lol)]- Ap- 
plying standard input-output formalism [l9l |. a°"*(ti;) = 
^yT\[cyS{uJ-iUL)+Sy{uJ)]-[c™S{uJ-uJL)+a'y], the output 
field is given by a°"'(^) = - "^;^'^" <(^) + iV^'-^v - 

d^)5{u~UL) + ^^^{iyiC,{uj-ujL)Hy2C2{u:-UL)), 
where dy is the input field of the probe laser. The back- 
ground terms Cy and dy only contribute for lo = lol. 
Homodyne measurements on the output field thus allow 



to measure ^yiCi(iu; — ujl) + £,y2C2{yJ — ^l)- The second 
probe laser on a mode with different and ^^2 measures 
a different linear combination of Ci and C2 and hence a 
different quadrature of the membrane vibrations. As the 
steady state of the membranes allows for repeated mea- 
surements, two probe lasers enable a reconstruction of 
the covariance matrix V, where the precision is limited 
by the input noise, a^". 

Conclusions - The scheme presented here allows to 
generate steady state entanglement of the motion of 
two dielectric membranes, which are suspended inside a 
Fabry-Perot cavity with a cavity decay rate that is lower 
than the mechanical resonance frequency of the mem- 
branes. The scheme only requires a drive by classical 
light and can work for environment temperatures up to a 
few Kelvin. With increasing environment temperatures 
stronger driving lasers and therefore shorter cavities with 
larger mode seperation are needed. The output fields of 
weak probe lasers that do not perturb the entanglement 
generation can be used to measure quadratures of the 
vibrational modes and thus enable a verification of the 
created steady state entanglement. 
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